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Abstract 

The complete calculation of the CMB polarization observables (i.e. E- and B-modes) is 
reported within the conventional ACDM paradigm supplemented by a stochastic magnetic 
field. Intriguing perspectives for present and forthcoming CMB polarization experiments are 
outlined. 



Large-scale magnetism recently became an intriguing triple point where cosmology, as- 
tronomy and high-energy astrophysics meet for complementary purposes [1]. Still we have 
no clues on its origin. The gravitational instability together with the subsequent galactic 
rotation could amplify a magnetic field spanning a collapsing region of the order of the Mpc 
in comoving units. The magnetic field regularized over such a scale L, i.e. Sl, is not em- 
pirically observable at the epoch of the gravitational collapse of the protogalaxy. But the 
Universe is a good conductor: the magnetic flux and helicity are approximately conserved 
implying that large-scale magnetic fields could have been already present at the time when 
photons last-scattered electrons and ions, i.e., according to the WMAP 5-year data [2], at a 
redshift z^cc — 1090. 

Intriguing effects related to tangled magnetic fields have been discussed with semi- 
analytical methods (see, in particular, [3]). More recently the impact of large-scale magnetic 
fields on scalar modes of the geometry have been addressed [4] and a dedicated numerical 
approach has been devised [5] . The complete calculation of the polarization angular power 
spectra (i.e., specifically, the EE, TE and BB angular power spectra) is here reported, for 
the first time, when the conventional ACDM paradigm is complemented by a stochastic 
magnetic field (i.e., according to the terminology of [5], mACDM scenario). 

In short the main theoretical impasse is the following. The large-scale description of 
temperature anisotropics demands a coarse grained (one-fiuid) approach for the electron- 
ion system: this is the so called baryon fiuid which is treated (with no exceptions) as a 
single fiuid in popular Boltzmann solvers such as COSMICS [6] and CMBFAST [7]. On 
the other hand the dispersive propagation of electromagnetic disturbances demands to treat 
separately electrons and ions, at least at high frequencies. It is appropriate to start from the 
Vlasov-Landau equation written in the form: 

^ + ^ 9^ ± + ^.^^^ )^ + 2^^.^ -Q^ - ^-11- 

where v = q/ \/im?a? + is the comoving three velocity, q is the comoving three- momentum 
and T is the conformal time arising, in the line element, as ds^ — a^(T){dT^ — [5ij — 
hij{x,r)]dx^dx^}. The prime denotes a derivation with respect to r. The rescaled elec- 
tromagnetic fields are denoted as E = a^E and as i? = o?B. By choosing the plus (minus) 
sign in Eq. (1), the evolution equation for the one-body distribution function f±{x,q,r) of 
the ions (electrons) can be obtained^ . In the electron-ion system Ccou is provided by Coulomb 
scattering. 

In the hmit e — > 0, Eq. (1) describes the evolution of neutral species. If Ccou = 0, Eq. 
(1) leads, below the MeV, to the well known evolution equation for the reduced phase space 



^The velocity-configuration space naturally arises since ions and electrons are all non-relativistic. Con- 
sequently, the comoving three- momentum is given by ^ = mav for each of the two charged species. The 
quasi-equilibrium distribution for electrons and protons is MaxwcUian and the strength of Coulomb scat- 
tering guarantees Tc ~ T, ~ T. For relativistic (neutral) species ~ rfq. The equilibrium distribution for 
neutrinos and photons will be, respectively, Fermi-Dirac and Bose-EInstein. 
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distribution of the neutrinos in the synchronous gauge^ [4, 5] : 



T'^ + ik\xT^ = 2//2(/i' + 6^') - 4^', n^k-n. (2) 

where, as in the conventional ACDM models the neutrinos are massless and, consequently, 
■y* = gYl^l = n*. The evolution equations of the brightness perturbations of the intensity 
(i.e. Ai) and of the polarization (i.e. Aq and Au), can be derived from Eq. (1) (always in 
the limit e — 0) by identifying Ccoii with the electron-photon collision term when the energy 
of the photons is parametrically smaller then the electron mass and when the electron recoil 
is neglected: 



(3) 



A'q + {tki, + 6')Aq = 1^(1 - l^')S, A^ + {iki, + e')Au = 0, (4) 

where e' — aXe "he crih is the differential optical depth defined in terms of the ionization 
fraction x^, of the electron density, and of the Thompson cross section (Txh; vi, — 9h/{ik) 
is the baryon velocity to be defined in a moment. The source term in Eqs. (3) and (4), 
i.e. S = (Ai2 + Aq2 + Aqo), contains the quadrupolc of the intensity of the radiation field, 
Ai2. In a nutshell, the CMB is linearly polarized (i.e. Aq ^ 0) since the amount of linear 
polarization is proportional, to first-order in the tight-coupling expansion, to the quadrupole 
of the intensity which is, in turn, proportional to the first-order dipole. This reasoning 
can be generalized to include the effects of the magnetohydrodynamical Lorentz force. The 
one-body distributions for electrons and ions enter Maxwell's equations as 

V • ^ = 47re y dMf+i^, V, r) - /_(f , v,t)], V • S = 0, (5) 
VxE + B'^0, V X B - E' ^4:ne J d\v[f+{x,v,T) - f.{x,v,T)]. (6) 

For length scales much larger than the Debye scale^, and for comoving frequencies smaller 
than the plasma frequency (i.e. U -C oJpe), Eqs. (2), (5) and (6) reduce to an effective 
one-ffuid description where the relevant dynamical variables are given by the centre of mass 
velocity of the electron-ion system (i.e. V}~, = {nicVc + 'm\v{)/{mc + mi)) and by the total 
current J. Ions (i.e. protons) are much heavier than electrons: the obtained equations can 
be expanded in powers of me/rrip <C 1. These approximations lead to (resistive) magneto- 
hydrodynamics (MHD) where the electromagnetic disturbances cannot propagate: the total 

— * — * — * 

current is solenoidal (i.e., according to Eq. (6), An J — VxB) and the electric fields van- 
ish in the baryon rest frame with an accuracy determined by the inverse of the Coulomb 



^Recall, for this purpose, that the scalar fluctuation of the geometry hij{x,T) carries two degrees of 
freedom, i.e., in Fourier space, hij{k,T) = [kikjh{k,T) + 2^{k,T){3kikj — %)]. 

^The Dcbyc length is defined as Ad = -\/T/(87re^no) with no = a^he = a^fn where fie and fn are the 
electrons and the ions concentrations. 
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Figure 1: The TE and EE angular power spectra. 



conductivity. The Lorentz force affects anyway the dynamics of the photon-baryon fluid: 

^7 = - 2^7b; '^^ ^ Y ~ ^ 4p^' 

where, by definition, 6*^ = V ■ t/x is the three divergence of the velocity field and 6^]^ = 6^ = 
6'b. In Eq. (8) 6^ and 5b are the density contrasts of photons and baryons. Both the metric 
fluctuations and the density contrasts of the various species will enter the corresponding 
Einstein equations whose explicit form can be found in [4, 5]. 

The comoving (angular) frequency corresponding to the maximum of the CMB spectrum 
is cJmax = 27rz/max wherc i/jnax = 222.617 GHz. The comoving plasma frequency is instead 
cJpe = 0.285 MHz for /i^fibo = 0.02273 (as implied by the best fit to the WMAP 5-year data 
alone). For CMB photon frequencies ZJ > ZUpe, Eq. (1) cannot be reduced to a one-fluid 
description. The stochastic magnetic field (obeying the one-fluid MHD equations) can then 
be treated as a background field. The propagation of the electromagnetic disturbances is 
calculated by taking into account the dynamics of ions and electrons within the cold plasma^ 
approximation which is rather safe since Tj -C mi and -C rrie and (^piasma ^ 1- The initial 
conditions of the Einstein-Boltzmann hierarchy are given in terms of the so-called magnetized 
adiabatic mode [4, 5] which is a solution of the system formed by Eqs. (7)-(8) (as well as by 
the other MHD and Einstein equations) in the tight-Thompson coupling approximation. The 
two sources of inhomogeneity of the system are represented by the stochastic magnetic field 
(which follows the effective set of one-fluid equations) and by the curvature perturbations. 



*The plasma parameter [8] is ^piasma = {Vuno)^^, i.e. the inverse of the number of charge carriers inside 
the Debye sphere. Around decoupHng gpiasma — 2.3 x 10"^^^ for the typical value of the baryonic density 
implied by WMAP 5-year data. 
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Figure 2: The BB angular power spectra induced by the dispersive propagation of the 
electromagnetic signal in a magnetized background. 

The Fourier amplitudes of the large-scale magnetic field will satisfy: 

{B,ik)B,ip)) = ^VBik)P,,ik)5^^\k + p), ^3 = ^3^-] . (9) 

where is the spectral amplitude, is the spectral index and /cl is the magnetic pivot 
scale; Pij = {6ij — kjij) the transverse projector. The curvature perturbations will be 
assigned consistently with the notations of Eq. (9) and, in particular, their power spectrum 
will be given by Vnik^) = •A-nik j kpY^~^ where An is the spectral amplitude at the pivot 
scale fcp = 0.002 Mpc"^; is the spectral index^. Our numerical code extends the code 
described in [5] and it is based on CMBFAST [7]. As decoupling approaches 0^ ^ 6'b and the 
linear polarization is generated. For frequencies of the CMB photons much larger than ciJpe, 
dispersive effect come into play (i.e. 9^ 7^ 9i) and the linear polarization is rotated^ leading, 
ultimately, to the BB angular power spectrum. Such a rotation is proportional ton- B where 
n, as before, is the direction of the photon momentum. The Larmor radius of the electrons 
and of the ions is much larger than the inhomogeneity scale of the magnetic field: the 
dynamics of electrons and ions (as well as the dispersion relations [9]) can be studied under 
the guiding centre approximation pioneered by Alfven [10]. The two helicities composing 
the (linear) CMB polarization propagate with different phase (as well as group) velocities. 
The Faraday rotation rate depends upon the difference uj[n^{uj) —n_{uj)]/2 where n±{lj) are 
the refractive indices of the two circularly polarized waves (one with positive helicity and 
the other with negative helicity). The heat transfer equations are then supplemented, in this 

^For the ACDM paradigm the 5-year WMAP data (alone) imply = 0.963t°:oi5. 

^There is a second class of dispersive effects which is related to the so-called ordinary and extraordinary 
waves [8]. The latter dispersion relations are insignificant since the scales of the problem imply that the 
refractive indices are 1 both for the ordinary and extraordinary waves [9] . 
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Figure 3: The bound on the BB angular power spectrum from the WMAP 5-year data (plot 
at the left). At the right, the magnetized EE and BB power spectra are compared to the 
B-modes expected from lensing and from the tensor mores in the case r = 0.3. 

regime, by the Faraday rotation rate: 

A'q + n'diAci = 2e'F(n)Au, A'y + n'diA^ = -2e'F{n)AQ, (10) 

where F(n) = 3/(16 '7r^e)n ■ B/u"^ and u denotes here the comoving frequency. Since the 
Faraday rate depends upon a stochastic field^ it will also be characterized by a power 
spectrum whose explicit form depends upon the spectral amplitude and slope of the mag- 
netic field [12, 9]. In Fig. 1 the EE and TE angular power spectra are reported in 
the case of the WMAP 5-year data. The TE and EE angular are defined from the cor- 
responding expansion coefficients by recalling that, in the present notations, A4±{n) = 
AQ(n) ±iAu(n) = a±2,^m ±2^m(^) where ±2^m(^) are the spin-2 spherical harmon- 
ics. In terms of a±2jrn the E-mode and the B-mode are given by = — (a2,^m + 0.-2, em) /"^ 
and by a)^ = i{a2,em — 0-2,£m)/2. In Fig. 1 the spectral index as well as all the other 
parameters of the underlying ACDM model have been fixed to the best-fit value of the 
WMAP-5year data alone. In both plots of Fig. 1 the best fit is illustrated with the full lines. 
The angular power spectrum of Faraday rotation is defined as C^^^ 6^1 mm' = {o*im^(-'m') 
where aim = I dQnY£rn{fi)*F{n). In terms of the autocorrelation of the B-mode will be 
given by 

= E ^2, ^)<^c^if ^, (11) 

^It is here assumed that spatial isotropy is unbroken (as observations seem to indicate). A uniform 
magnetic field (such as the one assumed in [11]) would break spatial isotropy. If the magnetic field breaks 
spatial isotropy the TB (and possibly EB) power spectra will be present. In the present case the latter power 
spectra vanish. 
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where Q{ii, £2, ^) is a function of the muhipoles containing a Clebsh-Gordon coefficient [12, 9], 
C^^^^ is the angular power spectrum of the polarization autocorrelations and^ 

CP - 30.03 He J^) "(1^)"- ^'^ + r (^M^ + ^ - i) 

V t^max / V ftL / p I "B-l j p / 6-nB j p / Z _|_ £ _ »B J 

where Qbl = -^l/ (S^rp^). In Fig. 2 the BB angular power spectra arc reported. In the plot at 
the left the dashed line shows the result obtainable from Eq. (11) in the case the Cf^^^ would 
be the one arising in the context of the ACDM adiabatic mode. The WMAP 5-year data (see, 
in particular, [2]) imply that, when averaged over £ = 2 - 6, £(£+ l)Cf^^^ / {2tt) < 0.15(/iK)2 
( 95% C.L.). The putative constraint of [2] does not make reference to a specific frequency. 
So it should be imposed at the lowest frequency channel. The lowest available frequency 
for this purpose would be for 27 GHz. The preceding frequency (i.e. 23 GHz) has been 
used as a foreground template and, consequently, the EE and BB power spectra have not 
been freed from the foreground contamination. We therefore choose to set the bound for 
a minimal frequency of 30 GHz since this is not only intermediate between the KKa and 
KQ bands of the WMAP experiment but it is also the putative (lowest) frequency of the 
Planck experiment [13]. In Fig. 3 (plot at the left) the full, dashed and dot dashed lines 
refer, respectively, to the cases oi u — 30 GHz, u — 100 GHz and u — Umax- In the right plot 
of Fig. 3 the magnetized EE and BB power spectra are compared with the B-modes from 
the lensing of CMB anisotropies and from the B-modes induced by the tensor modes (in the 
case of tensor-to-scalar ratio 0.3 which is the best fit value of the ACDM model plus tensors 
to the WMAP 5-year data). As it is apparent from Fig. 3 the constraints on the B-mode 
are not stringent for the magnetized background and are safely satisfied by a nG field at the 
epoch of the gravitational collapse ^. The constraints on the height of the acoustic peaks are 
comparatively more stringent [5, 9]. 

The obtained results suggest that multifrequency measurements of the GMB temperature 
and polarization within different channels will permit, for instance with Planck [13], an 
accurate scrutiny of the possible presence of magnetized birefringence. For this purpose, the 
scaling properties of the temperature and polarization autocorrelations in different frequency 
channels should be analyzed and compared. While should scale as , the EE and 

TT power spectra will be frequency independent [9]. 

K.E.K. is supported by the "Ramon y Cajal" program and by the grants FPA2005-04823, 
FIS2006-05319 and CSD2007-00042 of the Spanish Science Ministry. 

*For illustrative purposes we will limit our attention on the case ns > 1; in this situation Ab = 
(27r)"B-ir((nB - l)/2)Bl 

^The WMAP collaboration reports also limits on axion-induced birefringence. As discussed in [9] these 
bounds arc obtained by assuming that the birefringence is independent on the frequency of the incoming 
polarization (which is not true in the case of the magnetic field). Furthermore, the rate of axion-induced 
birefringence is fully homogeneous which is opposite to the case considered here. 
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